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Abstract 

The supergraph technique for calculations in super symmetric gauge theories 
where supersymmetry is broken in a "soft" way (without introducing quadratic 
divergencies) is reviewed. By introducing an external spurion field the set of 
Feynman rules is formulated and explicit connections between the UV countert- 
erms of a softly broken and rigid SUSY theories are found. It is shown that the 
renormalization constants of softly broken SUSY gauge theory also become exter- 
nal superfields depending on the spurion field. Their explicit form repeats that of 
the constants of a rigid theory with the redefinition of the couplings. The method 
allows us to reproduce all known results on the renormalization of soft couplings 
and masses in a softly broken theory. As an example the renormalization group 
functions for soft couplings and masses in the Minimal Supersymmetric Standard 
Model up to the three-loop level are calculated. 
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1 Introduction 



Supersymmetry, if it exists, must be broken. The gauge theory with softly broken 
super symmetry has been widely studied. To break supersymmetry without destroying 
the renormalization properties of SUSY theories, in particular the non-renormalization 
theorems and the cancellation of quadratic divergencies, one has to introduce the so- 
called soft terms They include the bilinear and trilinear scalar couplings and the 
mass terms for scalars and gauginos. 

A powerful method for studying SUSY theories which keeps supersymmetry man- 
ifest is the supergraph technique P, It is also applicable to softly broken SUSY 
models by using the "spurion" external superfields [|I|, ||, ||. As has been shown by 
Yamada with the help of the spurion method the calculation of the (3 functions of 
soft SUSY-breaking terms is a much simpler task than in the component approach. 
Using the spurion technique he has derived very efficient rules which allow him to cal- 
culate the soft term divergencies starting from those of a rigid theory in low orders of 
perturbation theory. 

However, these rules need some modification since, as mentioned by Yamada 
there might appear some problems in the calculation of vector vertices because they 
contain chiral derivatives acting on external lines. 

In the present paper we develop such a modification of the spurion technique in 
gauge theories which takes care of the above mentioned problem. We formulate the 
Feynman rules and show that the ultraviolet divergent parts of the Green functions of a 
softly broken SUSY gauge theory possess the factorization property being proportional 
to those of a rigid theory. 

In general, we consider a softly broken supersymmetric gauge theory as a rigid 
SUSY theory imbedded into the external space-time independent superfield, so that all 
the parameters as the couplings and masses become external superfields. This approach 
to a softly broken sypersymmetric theory allows us to use remarkable mathematical 
properties of iV = 1 SUSY theories such as non-renormalization theorems, etc. We 
show that the renormalization procedure in a softly broken SUSY gauge theory can 
be performed in exactly the same way as in a rigid theory with the renormalization 
constants being external superfields. They are related to the corresponding renormal- 
ization constants of a rigid theory by the coupling constants redefinition. This allows 
us to find explicit relations between the renormalizations of soft and rigid couplings. 

Throughout the paper we assume the existence of some gauge and SUSY invari- 
ant regularization. Though it is some problem by itself, we do not consider it here. 
Provided the rigid theory is well defined, we consider the modifications which appear 
due to the presence of soft SUSY breaking terms. To be more precise, in the following 
sections when discussing one, two and three loop calculations of the renormalization 
constants we have in mind dimensional reduction and the minimal subtraction scheme. 
Though dimensional reduction is not self-consistent in general, it is safe to use it in low 
orders and all the actual calculations are performed in the framework of dimensional 
reduction H, ^ Nevertheless, our main formulae have general validity provided 
the invariant procedure exists. 
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2 Softly Broken Pure N = 1 SUSY Yang-Mills 
Theory 



Let us first consider pure = 1 SUSY Yang-Mills theory with a simple gauge group. 
The Lagrangian of a rigid theory is given by 

C„gid = J d^e -^TiW'^Wa + J d^e -l^TrW^Wa. (1) 

Here W"' is the field strength chiral superfield defined by 

where = (-^'^)^ ^'^'^ generators of the gauge group G. 

To perform a soft SUSY breaking, one can introduce a gaugino mass term. This is 
the only term allowed in a pure gauge theory which does not break a gauge invariance. 
The soft SUSY breaking term is 

-r -H^w^II^u 

so ft— breaking 2 2 ' 

where A is the gaugino field. To rewrite it in terms of superfields, let us introduce an 
external spurion superfield 1] = 6^, where ^ is a grassmannian parameter. The softly 
broken Lagrangian can now be written as 

Csoft = J d^e ^(1 - 2iie^)TiWw^ + j d^e ^(i - 2iJie^)TiWw^. (2) 

In terms of component fields the interaction with external spurion superfield leads to a 
gaugino mass equal to itla = /U, while the gauge field remains massless. This external 
chiral superfield can be considered as a vacuum expectation value of a dilaton superfield 
emerging from supergravity, however, this is not relevant to further consideration. 



2.1 The Feynman rules 

Consider now the Feynman rules corresponding to the Lagrangian (0). For this purpose 
it is useful to rewrite the integral over the chiral superspace in (j^) in terms of an integral 
over the whole superspace. One has 

Csoft = j d^Od^e (^{1 - 2/i^2^Tre-^D"e^r)2e-^D,e^ 
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+^(1 - 2fie^)TTe-^D''e^D^e-^D^e^j . (3) 



Consider the quadratic part of the action 

^(2) = J d^xd^e d^OTi ^ (l - D''V{x, 9, 9)D^D^V{x, 9, i 
+ J d^xd'^9 d^Ti ^ (l - 2/i^2^ D^V{x, 9, 9)D^D^V{x, 9, i 
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To find the propagator, one first has to fix the gauge. In a rigid theory the gauge 
fixing condition is usually taken as 



D'V = f, D'V = f. 



(4) 



This is equivalent to adding to the action (|l]) a gauge-noninvariant term proportional 
to 



1 ^ 
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■■gauge— fixing Tr ^ f f Tr 



(5) 



By joining this together the quadratic part of the rigid action takes the form 



am 

rigid 



It can also be rewritten as 



9' 



ni/2 + -Ho 



a 



where we have introduced the projectors 



ni/2 



8^2 ' " 16^2 
The inverse operator then gives the vector propagator 

_ ni/2 + aUo 



V^iz,)V\z2) 



rigid 



Q2 



ah 



The next task related to the determination of the vector propagator is to find the 
associated ghost term. Under the gauge transformation the superfield V is transformed 
as 



or in the infinitesimal form 



V -^V + H{V)A + H{V)A 

^^A[(A + A)+cothQA(A-A 



V 



Consequently, the associated ghost term is 



Cghost = J d^O Tr ^bD^ {H{y)c + H{y)c) + j d^ Tr ^6^^ ^h{V)c + H{V)c 



/ciWTrl(6 + 6)(4)A 



[c + c) + coth(^) A(c-c) 



(6) 



where b and c are the Faddeev-Popov ghost chiral superfields. 
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For a system in the external spurion field we can change the gauge fixing condition 
(Q) in order to take into account the interaction of the associated ghost fields with 
the spurion. The reason for this is the appearance of these terms in due course of 
renormalization. Here, we depart from the Feynman rules of ref.Q. Our aim is also to 
get the same common spurion factor for the ghost propagator as for the vector one, as 
it will be clear below. For these purposes we choose the gauge-fixing condition in the 
form 

D'-Al- l^Q' - l^Q'^ V = f, D'^(l- ^,9' - ^,6') V = f. 

or 



g2 \ J gl 



9 9 
where we have introduced the notation 



D^-V = f, D'-V = f. (7) 



= g^(l + ^^e^ + + 2flW) . (8) 
Note that due to the grassmannian origin of 6 

[i - fio^ - fxe^) (i + fi9^ + fiO'^ + 2fx^e^e^) = i. 

Now the gauge-fixing term can be written in the standard form 



gauge— fixing 



I dVOTi - ff, (9) 



r 

where ^ might be a real constant superfield. In what follows we take 

^ ~9'' 

which leads to the following quadratic part of the action in a softly broken theory 



^ 'l~2iie'')D''VD''D^V (10) 



One now has to find the inverse operator. It can be done analogously to the rigid 
case by introducing the projection operators. One should take into account, however, 
that while integrating the covariant derivatives by part they may also act on the spurion 
field. Strictly speaking, this leads to additional terms in the propagator which are, 
however, suppressed by powers of momenta and are inessential in the analysis of UV 
divergences. To understand it better, one can treat the spurion terms in eq.(p!OD as 
interactions, and to find the propagator, one has to take into account an infinite chain 
of this kind of insertions adding up into a geometrical progression. It can be done by the 
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method of Ref . . If in due course of this procedure some of the covariant derivatives 
act on a spurion field, one does not get enough powers of momenta to cancel the 
denominator and the resulting terms are thus suppressed for a high momentum square. 
Thus, as it has been also noted in Ref. P] , for the calculation of the divergent parts of 
the Green functions it is sufficient to take the part of the vector propagator where the 
chiral derivatives do not act on spurious. 

Having this in mind, one can proceed in full analogy with the rigid theory, replacing 
the coupling by and get 

V\zi)V\z2)) = g! ni/2 + «no ^(8)^^^ _ ^ irrelevant terms, 

/ soft 4 (j^ 

where by irrelevant terms we mean the ones decreasing faster than for large p^. 
Thus, we get a simple relation between the vector propagators of soft and rigid theories: 



F(xi,^^l,^^l)\/(X2,^2,^2)) , = (11) 

/ soft 

= ^ (V(xi,6i,9i)V (x2, 02,02)) + irrelevant terms, 

^2 \ / rigid 

Notice that contrary to ref. [^] this relation is valid for any choice of a. 

The change of the gauge fixing condition leads to the change of the ghost Lagrangian 
(|^). It becomes 

CgHost = Jd''OTTbD^(^^{H{V)c + H{V)c)j 

+ J d^O Tr ID^ {H{y)c + H{y)c)^ 

= /rfWTrl(6 + 5)(-^A[(c + ^)+cothQA(^-s)])- (12) 

As one can see from eq. (0), the situation with the ghost propagator and the ver- 
tices that contain ghost superfields is more simple. This is due to the absence of chiral 
derivatives in the ghost-vector interaction. Nevertheless, when calculating the inverse 
operator one has to repeat the same procedure with the covariant derivatives arising 
from the source terms. In complete analogy with the case of the vector propagator 
the essential part of the ghost propagator is obtained by assuming that the covariant 
derivatives do not act on the spurion field. The other terms are suppressed by the 
powers of momenta. Hence, one again has 

(G(zi)G(z2)) = (flg^) Ig{zi)G{z2)) + irrelevant terms, (13) 

\ / soft \ / \ / rigid 

where G stands for any ghost superfield. 

Thus, to perform the analysis of the divergent part of the diagrams in a soft theory, 
one has to use the same propagators as in a rigid theory multiplied by the factor j . 
It is also obviously true for any vertex of the ghost-vector interactions of the softly 
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broken theory (|T2]). Each vertex of this type has to be muhiphed by the inverse factor 

The situation is less obvious with the vector vertices. Here, one has two terms 
as in eq.(0) which differ only in the order of covariant derivatives, one being complex 
conjugated to the other. Performing grassmannian algebra in any diagram, one can 
always replace one term by the other; so actually one has only one term. This is 
what we have in a rigid theory. In a soft theory the situation is similar but the first 
term is multiplied by a factor (1 — , while the second by the complex conjugated 
one 

(l-2^^2y Hence, one can consider the vector vertices of only one type, as in a 

rigid theory, multiplied by a factor (l — fiO'^ — jj,6'^^, i.e. by g'^/g'^ as in the case of the 
ghost-vector vertices. 

Thus, we see that any element of the Feynman rules for a softly broken theory 
coincides with the corresponding element of a rigid theory multiplied by a common 
factor which is a polynomial in the grassmann coordinates. 

2.2 The ultraviolet counterterms 

Consider now the Green functions of a softly broken theory constructed with the help 
of the Feynman rules derived in the previous section. Due to the presence of a factor 
g^/g"^ in the vertices and of the inverse one in the propagators, any Feynman diagram 
in a broken theory is given by that of a rigid theory with the spurion factors on the 
lines and vertices. Our aim is to show that these polynomials in grassmann variables 
can freely flow through the diagram and are finally collected in front of a diagram as 
a common factor. 

Indeed, consider the corresponding grassmannian integration. The usual procedure 
of evaluation of grassmannian integrals in a rigid theory is the following: 

i) If there are no chiral derivatives Da or Da on the line, it is proportional to the 
(5-function and can be shrunk to a point. 

ii) If there is some number of derivatives, one can remove them from a line integrat- 
ing by parts. 

iii) Repeating this procedure several times and using the properties of chiral deriva- 
tives 

{Da, Df,} = {Da, Dfs} = 0, DaDf, = eapD\ D^D^D^ = IGp^D^, etc. 

to reduce their number, one can finally come to a final integration where all the 
chiral derivatives are concentrated on one line. The powers of momenta which 
appear in due course of this procedure cancel some propagator lines so that the 
resulting diagram diverges. 

iv) The final integration is performed with the help of the relation 

6uD^D^u = 165i2. (14) 
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Consider what happens if one adds grassmannian polynomials on the lines and 
vertices. Suppose first that all chiral derivatives inside a graph are contracted and 
integrated inside the graph so that there are no external derivatives left. If one of these 
internal derivatives acts on a spurion field it creates additional 9°" instead of Z)". Then 
either this 9°" survives till the last integration, or some other Da acts on it. Thus, at 
the last step of integration one has instead of either 

S^29^DlD'^5u or 6i2Dl5i2. 

It is easy to see that both the expressions are zero. Thus, the configuration when the 
chiral derivative acts on a spurion field gives no contribution. This means that the 
spurion field flows through the diagram freely and can be factorized. 

Consider now the case when a graph contains some derivatives that act on external 
lines when performing the grassmannian integration. This case corresponds to vector 
interactions. 

Let some covariant derivative, which originally was inside the diagram, pass through 
the diagram during the above mentioned procedure and act on the external vector line. 
We can trace this derivative in the diagram of the soft theory. On its way through the 
diagram it can meet some "spurion" factor in a vertex or a propagator. Integrating by 
parts one has two terms. The first one, when the chiral derivative acts on a spurion, 
and the second one, when it acts on the neighbouring propagator. When the covariant 
derivative does not act on a spurion it goes through the diagram like in a rigid theory 
and acts on external line. As a result one gets the expression {DaV"^). In the 
other case, when the covariant derivative acts on the spurion, the number of covariant 
derivatives is reduced compared to the rigid theory. Then, either one does not have 
enough of them to cancel the powers of momenta in the denominator or to get the 
non-zero answer or one has the same diagram as in the rigid case, but with the factor 
Dg^, that goes outside the diagram. Independently of the position of this factor inside 
the diagram, the final expression one gets after the calculation of the diagram contains 
[Dag^) V"^ which together with the usual contribution, gives Da {g'^V"^) instead of 
g^DaV"^ in the rigid case. 

All together these factors g"^, which are met by the chiral derivative on its way to 
the external line of the diagram, are collected in a monom R{g'^). Therefore, the final 
expression for the diagram contains Da {R{g^)V'^) . The same conclusion is true for all 
four derivatives of the divergent vector diagram that act on the external lines. Thus, 
the counterterm to this diagram of the soft theory is 

R,{f)V^D" (R2{f)V-) D' [R,{f)V^DaR,{f)V') , (15) 

where i?2, -Rs and R^ are some monoms in g"^. This expression should be compared 
with the counterterm for the corresponding diagram of the rigid theory 

R{g^)V'^D''V''D^V''DaV\ 

where 

R,{g')R2{g')R,{g')R,{g^) = R{g'). 
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The momentum integral for the diagram in a soft theory is identical to that in a rigid 
one since the algebraic operations with the covariant derivatives coincide for the both 
cases. 

We now want to argue that the monoms Ri{g'^) can stand anywhere in eq. (|l5|) , 
i.e. that the spurious can flow through the diagram. Indeed, in a rigid theory the 
counterterms are absorbed into the redefinition of the field V and the coupling g"^ 

Vb = Zl/\g')V, 
gl = Z,ig')g' 

In a soft theory Ri are the functions of the external spurion superfield and the renormal- 
ization constants become these functions too. Collecting all counterterms of the type 
([l5|) , one should get an expression which being written in terms of the bare quantities 
obeys the gauge invariance, namely 

Tr R[{f)e-''-D'^RXg')e''^D' (R',{f)e-''-D^R',{f)e''-) , (16) 

where the primed monoms mean that partly they are absorbed into the renormalization 
of the field Vb. 

Now it is clear that the covariant derivatives in eq.(|T6D do not act on monoms 
Ri{g'^). Indeed, if one of the covariant derivatives in eq.(|IBp acts on the momons 
R'2 or i?4, the two exponents cancel each other and one is left with Tr that is equal 
to zero. In the case when the both chiral derivatives in eq. (|TB|) act on R'2 and R'^ 
in both cases the exponents cancel and one has a divergent constant that should be 
absorbed into the renormalization of the vacuum density. 

Hence, one can take the monoms R'2 and R'^ out of the covariant derivatives Da. The 
other derivative, i)^, does not act on the monom for the reason that e~^^ = 
so that in the bracket in eq.([T6|) only the chiral spurion rj = 6"^ survives and does 
not act on it. 

Thus, we come to a conclusion that the monoms can be factorised in front of the 
expression in eq. (|T^) just like in the rigid theory and the only difference between the 
rigid and softly broken theories is that the coupling constant g^ should be replaced by 
g^, i.e. 

Zi = Zi [g'^ , 

3 Softly Broken SUSY Gauge Theory with Chiral 
Matter 

Consider now a rigid SUSY gauge theory with chiral matter. The Lagrangian written 
in terms of superfields looks like 

Crigid = J d^Od^e ¥{e^)l<i>, + Jd^e W + J d^9 W, (17) 
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where the superpotential W in a general form isQ 



W = lx'''<^,<^j<^k + ^M^^$,$,-. (18) 
6 2 



The SUSY breaking terms which satisfy the requirement of "softness" can be written 

as 

^soft— breaking 



+ (m2);.0*0^ (19) 



Like in the case of a pure gauge theory the soft terms (|TUp can be written down in 
terms of superfields by using the external spurion field. Thus, the full Lagrangian for 
the softly broken theory can be written as 

Csoft = JdW ¥i6^-im')1vv)ie'')i<!>, (20) 

+ I (fe ^(A'^'^ - A'^''r^)^,^,^k + \{M'^ - B'^ri)^i^j + h.c. 

3.1 The Feynman rules 

The Lagrangian (^) allows one to write down the Feynman rules for the matter field 
propagators and vertices in a soft theory. We start with the propagator of the chiral 
field. For the purpose of the analysis of divergences one can ignore the mass terms M*-' 
and 5*-^, since in the minimal scheme they do not contribute to the UV divergences. 
Then, the quadratic part looks like 

The inverse operator is easy to obtain due to the nilpotent character of the spurion 
fields. One easily gets 

{^{z,)Mz2y) = (5f + ($ (.i).«i>(^2)')„^^^ (5/ + \{m')\rin) 

+ irrelevant terms, (21) 

where the irrelevant terms again arise from covariant derivatives acting on spurious 
and are suppressed by powers of momenta. 

The vector- matter vertices, according to eq. (pO]) , gain the factor {5l — {m'^)lrif]) so 
that if in a diagram one has an equal number of chiral propagators and vector-matter 
vertices the spurion factors cancel. 

The chiral vertices of a soft theory, as it follows from eq . (pOf ) , are the same as in a 
rigid theory with the Yukawa couplings being replaced by 



^ To avoid further complications (see, e.g. ||^) we do not introduce linear terms in eqs.(^ 19). 
This corresponds to the case when none of the fields goes into vacuum. In the absence of singlets this 
is guaranteed by the gauge invariance imposing some restrictions on the couplings. 
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3.2 The ultraviolet counterterms 



The structure of the UV counterterms in chiral vertices is similar to that of the vector 
vertices, but is simpler due to the absence of the covariant derivatives on external lines. 
This corresponds to the first case considered in the previous section. To get the UV 
divergent diagram, the covariant derivatives should not act on spurion fields, which 
means that spurious factorize. In the diagrams, when the chiral vertices are contracted 
with the chiral propagators, this results in the following effective change of the Yukawa 
couplings of a softly broken theory 

^ijk _^ yjk ^ yijk _ j^ijk^ ^ }_i^ynjk^^2-^^ ^ A*"'=(m2)^„ + X^''{m^)l)7]f], (22) 

Kjk Kjk = Kjk - AjkV + ^{Kjk{m^)7 + \nk{m^Tj + Kjn{m^)k)vV- (23) 

Thus, the UV counterterms of a softly broken theory are obtained from those of a 
rigid one by the substitution 

g"^ g^, ^^^^ \jk Kjk- 

4 Renormalization of Soft versus Rigid Theory: 
the General Case 

The external field construction described above allows one to write down the renormal- 
ization of soft terms starting from the known renormalization of a rigid theory without 
any new diagram calculation. The following statement is valid: 

The Statement Let a rigid theory (^\T3J be renormalized via introduction of 
the renormalization constants Zi, defined within some minimal subtraction massless 
scheme. Then, a softly broken theory ^j^) is renormalized via introduction of the 
renormalization superfields Zi which are related to Zi by the coupling constants redefi- 
nition 

Ug,X,\...) = Z,if,~X,\...), (24) 
where the redefined couplings are 

f = g'^i'^ + i^v + I^v + '^mv), ^ = 6'^ ^ = ^^ (25) 

hjk = Kjk - Ajkf] + ^{Xnjkirn'^Yl + Knkim'^yf + Kjnim^^ (27) 

This allows us to find explicit relations between the renormalizations of soft and rigid 
couplings which is an explicit realization on the level of final expressions of the rules 
A and B of ref.|]. 

The renormalization constants for vector superfields and gauge couplings are general 
superfields, i.e. they depend on t] and f/, while those for the chiral matter and ghost 
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fields and the parameters of a superpotential are cliiral superfields, i.e. tliey depend 
eitlier on 77 or f/. 

From eqs.p^) and (p5| - p^ ) it is possible to write down an explicit differential oper- 
ator which has to be applied to the j3 functions of a rigid theory in order to get those 
for the soft terms. We first construct this operator in a general case and compare the 
resulting expressions with explicit calculations made up to two loops. Then, using the 
formulated algorithm we calculate some three-loop soft term j3 functions. In the last 
section, we consider some particular models. 

To simplify the formulas hereafter we use the following notation: 



4.1 Pure gauge theory 

Consider first the gauge couplings a^. One has 

^Bare rv ^ -Bare ry ~ ( nQ\ 

where Z^i is the product of the wave function and vertex renormalization constants. 

Though <5j and Z^i are general superfields, one has to consider only their chiral 
parts since the Lagrangian consists of two terms, a chiral and an antichiral, and in 
each term only the proper chirality part contributes. Therefore, we consider the chiral 
part of eq. (|28|) 

+ mj'-^r;) = + mA,7])Z^,{a)\^=o- 

Expanding over rj one has 



«f"^^ = (29) 

^Bare^Bare ^ m A^a^Z ^,{a) + a.D^Z (30) 



where the operator Di extracts the linear w.r.t. r] part of Z^i^a). Due to eqs. (p5| - p^ ) 
the explicit form of Di is 

Di = rriAiOii— . 

dai 

Combining eqs.(p9|) and (|30| ) one gets 

m'i^'' = mA,+D^\nZ^i . (31) 
To find the corresponding j3 functions one has to differentiate eqs.(^) and (pl|) 



w.r.t. the scale factor having in mind that the operator Di is scale invariant. This 
gives 

where jai is the logarithmic derivative of In Z^i equal to the anomalous dimension of 
the vector superfield in some particular gauges. 



This result is in complete correspondence with that of Ref. |14|. Indeed, from 



equations (52) one can derive that the ratio mACt/Pa is renormalization group invariant. 
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4.2 Chiral matter 

Consider now the chiral matter. Due to the non-renormahzation theorems, there is 
no renormahzation of a superpotential. This means that the only renormalization 
comes from the wave functions. The corresponding term in the Lagrangian including 
renormalizations looks like 

" (fed^e ¥z{^j, (33) 
where the renormalization superfield now has a decomposition 

Z]=~z\^(b\^^\m)z\,. (34) 

where the chiral(antichiral) renormalization superfields Z\^^{Z^^ are the wave function 
renormalizations and A* is the soft term (m^)* renormalization. To find, them expand 
Zj{a,'y,y) over the grassmann variables 

Zj = Z] + D^Z]!] + D^Z]fi + ^2^;^, (35) 

where the operator Dx now has to take into account the ri dependence of the Yukawa 
coupling y"^^^ I)\ is conjugated to D\ and is a second order differential operator 
which extracts the 777/ dependence of Zj. They have the form 

= (37) 

h-nr brv 9 8 8 8 \ 



Qybac Qybca Qy^^^ Qy^^^ 8ybcn) 

Substituting this into eq. (|3^) one has 

Z;^. = Zl^ + {Z^%D,Z^i^ = Zl,{5'; + {Z-'f,D,Z]r^), (3J 



z^^ = z;^ + D,zi{z^'f^n = {^l + D,zi{z-')[n)z'^^, (39) 

A} = {Z~^%DX{Z-^')] - {Z-,%D,Z'^{Z-XDiZI{Z-,X. (40) 
where Z^^- satisfies 

and is the square root of in the perturbative sense. The inverse one is 

{Z^% = {61 - {Z-^)\D,Zir^){Z^^)]. (41) 

We can now write down the relations between the renormalized and the bare cou- 
plings of a superpotential. One has 

{MBareT = {Z^%{Z^yi , (42) 

{ySarer' = y'^"- {Z^')\{Z^')UZ^'t (43) 
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or 



and analogous for y^^''. Expanding over t] one has for M*-' and B^^ 
M'Le = M^\Z~,%{Z-^^)i, 

Btre = B'^\Z-,%{Z-,')l 

+M^\Z^%{Z-y^D,Zl^{Z^^)T + M'\Z-^)\^D,Z^{Z^')l{Z^')\. 
The last line can also be rewritten as 

= B^\Z-/)l{Z^^)l + M|_(Z-iDiZ)i + Ml^_(Z-iDiZ)| (44) 
Differentiating with respect to a scale factor one gets after some algebra 

^ = liMSi + M^^Yi), (45) 
P'i = \{BSi + B'^^]) - {M^'D.^l + M'W,^\), (46) 



where 7* is the anomalous dimension of the matter chiral field being the logarithmic 
derivative of ln(Z~^)*. 

Analogously, for the Yukawa coupling y^^^ and for the soft triple coupling A^^^ we 
obtain 

/3f = \{y'H + y'''il + y'''i]). (47) 

/3f = \{A^H + A'^'ll + A'^^'li) (48) 
Consider now the renormalization of the soft mass term {m'^)'j. One has 

«are)] = {m')] " A} (49) 

with A* given by eg. (PH) . Differentiating it with respect to a scale factor and using 
the explicit form of the second order differential operator D2, after some algebra one 
obtains a simple expression for the (3 function 

(AnO} = ^27} (50) 

4.3 Summary on the soft term renormalizations 

Relations between the rigid and soft terms renormalizations are summarized in the 
Table. 
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The Rigid Terms 


The Soft Terms 






D2 = D,Di + m%ai^^ 

1 l/'™2^a I „,nbc 9 , „,bnc , r,.bcn , „, 9 , „, , „, 9 \ 



Table 1: Relations between the rigid and soft term renormalizations in a massless 
minimal subtraction scheme 

5 Illustration 

To make the above formulae more clear and to demonstrate how they work in practice, 
we consider the renormalization group functions in a general theory up to two loops. 
We follow the notation of ref . [§] except that our j3 functions are half of those of ref . [§] . 
Note that all the calculations in ref. p| are performed in the framework of dimensional 
reduction and the MS scheme. 

The gauge (3 functions and the anomalous dimensions of matter superfields in a 
massless scheme are the functions of dimensionless gauge and Yukawa couplings of a 
rigid theory. 

5.1 One- loop renormalization 

In the one-loop order, the renormalization group functions of a rigid theory are (for 
simplicity, we consider the case of a single gauge coupling) 

7i^) = «Q, Q = T{R)-3C{G), (51) 

7-^'^ = ly'%M-2aC{R)], (52) 

where the Casimir operators are defined by 

T{R)6AB = Tr{RARB), C{G)6ab = lACDfBCD, C{R)) = {RaRa)). 
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Using eqs.(|32|-^), we construct the renormalization group functions for the soft 
terms 



= arriAQ, (53) 
1 



+M'i(^^A^krny^km + 2am^C(i?)/) + j), (54) 
2 ^2' 

^yVi^^A'^muy^^^ + 2«m^C(i?)f ) 

+ {i ^ j) + (i ^ k), (55) 
1 

2" 

+ iy"'=V')n%fc^ + ^y*'V)"l/nM + ^y''\m')>;y,ki. (56) 

One can easily see that the resulting formulae coincide with those of ref. ^ with one 
exception: we have ignored hereafter all the tadpole terms assuming that they all are 
equal to zero because of the absence of linear terms in the action (see the footnote on 
p.8). 

5.2 Two-loop renormalization 

In two loops the situation is more tricky. The rigid renormalizations are 

7i') = 2a2^(G)g-— C(i?)5(VWi-2aC(i?)|), r = dimG = Saa, (57) 

r ■'2 

^.(2) ^ -(y^^Py^^^ + 2aCiR)X)i^y-''%ki-2aCiR);) + 2a'QCiR)]^ (58) 

Then, the soft renormalizations are as follows: 

/3S = Aa'mAC{G)Q-^^CiR)]i^y^'^^y,,i-2aC{R)i) 

+—C{R)]{lA^'''yM + 2amAC{R)i), (59) 

PP''> = -\B^'{y''^yikn + 2aC{R)f6i){^y^^%,t-2aCiR);) 

-M^\A^^^yi,,, - 2amAC{RYi5i){Uj'^^%st - 2aC{R);) 

-M'^y^'^Pyikn + 2aC{R)^6i){^A'''%st + 2amAC{R);) 
+B'^a^QC{R)i - AM'^a^QC{R)\mA + {i ^ j), (60) 
Pf^'^ = -^A^^'(rt™„ + 2aC'(i?)f5^)(irtp.*-2aC(i?);) 
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+A'^^a^QC{R)'[ - Ay'^^a^QC{R){mA 

+ ]) + {i^k), (61) 

2 11 

+ _/p(rn2)^^^.^^^ + i^y"^"{m^)lyjkn + ^Vikpim^Ynyjks 

-{y'^'y^kn + 2aC(i?)X)(^^"^*^p.t + \{m')^y''%st 
+ ^^y-^^yUm% + ^y"'*(m2)fy,,, - 2am\C{Rr,) 
-{A^'^^y.kn - 2amAC{R)%){]^y''''A,,t + 2am^C(i?)^) 

-(y*'^^'A,fc„, - 2amAC(i?)^(5;)(^A"^Vpst + 2amAC{R);) 

+l2a^m\QC{R)] (62) 

Comparing these formulae with those of ref. |§] we find that they coincide with one ex- 
ception. The renormahzation of m? in eq.(16) of ref. p[ has two extra terms, 8g^SC{Ryj 
and the term proportional to the mass of the so-called e-scalars, m^. The first term, 
though it is not proportional to m^, still has the same origin. It has appeared as a 
result of the e-scalar mass counterterm in one loop. 

The presence or absence of these terms is renormahzation scheme dependent. The 
version of dimensional reduction adopted in ref. corresponds to diagram by diagram 
minimal subtraction of divergences and naturally includes the e-scalar mass countert- 
erm. Our substitution being formulated in the superfield formalism does not contain 
5-scalars; no surprise that the corresponding terms do not appear. 

There are two ways to resolve the noticed discrepancy. Either to redefine the 
renormahzation scheme in such a way that these terms are not present, like the one 
in Ref. |T^, or to modify the superfield renormahzation procedure to reproduce these 
terms (see also the discusson of this point in Ref. |[T5[| ). 

5.3 Three-loop soft renormalizations 

To demonstrate the power of the proposed algorithm, we calculate the three loop 
gaugino mass renormahzation out of a gauge (3 function. One has in three loops ^ 

7(3) = a^C{G)Q[4C{G)-Q]-^a^QC{R)iC{R)i 

+ -a\y'%M - 4aC{R)])C{R)iC{R)t - -a'CiG)iy''%M - AaC{R)])C{R)l 
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+ ^^aiy^'^'y.ki - AaC{R)i){y^^%,t - AaC{R)l)C{R)l (63) 
The corresponding gaugino mass renormahzation is 

/^S = ^a^C{G)Q[AC{G) - Q]mA - ^a^QC{R))C{R){mA 

+ ^a\f%M - 4aC{R)))C{R)iC{R)tmA 

- ^a'iA'%M + AaC{R))mA)C{R){C{R)l 

- ^a'C{G){y'^'y,ui - AaC{R)])C{R)lmA 
+ ^-a'C{G){A'^%ki + ^aG{R))mA)C{R)l 

+ ^^ay^'''^y,Uy''''ymst - AaG{RrjG{R)lmA 

- ^^aA'^'^y.Uy^^'ymst - AaG{RrjG{R)l 

- |:a2/''"%fen(A"^*y^,t + AaG{R)l,mA)G{R){ 

+ ^^aiy'^'ym - ^aG{R)]){y^'%st - AaG{R)l)G{R)lmA 



^-a^A'^^'y^ki + AaG{R))mA){y'''%st - 4aG{R%)G{R)l 



- -a{y'%M - AaG{R)]){A^'%,t + 4aG{R%mA)G{R)l (64) 

6 Soft Renormalizations in the MSSM 

The general rules described in the previous section can be applied to any model, in 
particular to the MSSM. In the case when the field content and the Yukawa interactions 
are fixed, it is more useful to deal with numerical rather than with tensor couplings. 
Rewriting the superpotential ( [TE| ) and the soft terms (|T^ in terms of group invariants, 
one has 

^susY = ^ E 2/aAf + i E M,h^^,^,, (65) 

and 

— ^soft 



Q^.^-a .... 2 , 



^ E •Aa\a''(f)i(j)j(j)k + ^ E Bbhh(pi(pj + ^rUAjXjXj + h.c. 
^ b 



(66) 

where we have introduced numerical couplings ya, Mb, Aa and Bb- 



17 



Usually, it is assumed that the soft terms obey the universality hypothesis, i.e. they 
repeat the structure of a superpotential, namely 

Aa = yaAa, = M^B,, {m')]=ml6]. (67) 

Thus, we have the following set of couplings and soft parameters: 

Qj, Ua, Mb, Aa, Bb, m-, niAj- 

Then, the renormalization group (3 functions of a rigid theory (^,^,^) look like (for 
simplicity, we assume the diagonal renormalization of matter superfields) 

Paj = Pj = aj'Jaj, (68) 

Pya = \yaY.^^i^i^ (69) 
Pm, = iMbY^Tb^i, (70) 

where ji is the anomalous dimension of the superfield 7q,^ is the anomalous dimen- 
sion of the gauge superfield (in some gauges) and numerical matrices K and T specify 
which particular fields contribute to a given term in eg. (|65|) . 

To get the renormalization of the soft terms, one has to apply the algorithm of the 
previous section, eqs.( |5^ , ^ , |i5| ). In terms of numerical couplings it is simplified. 



The renormalizations of the soft terms are expressed through those of a rigid theory 
in the following way: 

PrriA. = Di-faj, (71) 

Pa. = -D^Y.Kaa^, (72) 

i 

Pb, = -D.Y^Tbau (73) 

i 

Pm^ = ^27., (74) 

and the operators Di and D2 now take the form 

A = -^-"^^ - ^'^^'^ J;' (75) 
A = (-^^»«^^-^'^^«Jr)^ + <«.A+^2^^^y„J^. (76) 

2 

where we have used the notation = 
6.1 Illustration 

To illustrate these rules, we consider as an example one loop renormalization of the 
MSSM couplings. Consider for simplicity the third generation Yukawa couplings only. 
Then, superpotential is 

y\^MSSM = {ytQ'U^m + VbQW'Hl + VrUE^Hl + fiHlHi)ei„ (77) 
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where Q,U,D', L and E are quark doublet, up-quark, down-quark, lepton doublet and 
lepton singlet superficlds, respectively, and Hi and H2 are Higgs doublet superfields. i 

and j are the SU (2) indices. 

The soft terms have a universal form 

— ^soft-breaking — X] "^i + ( o X] '^"'^a C''^) 

i ^ a 

+ Atytq^u'^h^ + Abybq^d"'h[ + ArVrPe'^h^ + Bi^h\h{ + h.c), 

where the small letters denote the scalar components of the corresponding superfields 
and Xa are the gauginos. The SU (2) indices are suppressed. 

Renormalizations in a rigid theory in the one loop order are given by the formulae 

33 

= biai, 1,2,3, bi^—,1,-3, (79) 

= Yt + Yb-^a3-^a2-^ai, (80) 

7?^ = 2Yt-^as-^ai, (81) 

7!,') = 2n-^a3-^ai, (82) 

7i'^ = Yr-^a^-^ai, (83) 

7^'^ = 2Yr--ai, (84) 
o 

7^ = 3n + n-^«2-^«i, (85) 

tS = 3y,-^a2-^«i. (86) 

Consequently, the renormalization group (3 functions are 

/32^ = kal (87) 

= y^(6y, + y,- Q;3-3a2--ai), (88) 
3 15 

16 7 

= n(y* + 6n + n-y«3-3a2-— ai), (89) 

= i;(3n + 4n-3a2-^ai), (90) 

o 

tf.^ = /x'(3Yt + 3y6 + y,-3a2-^ai). (91) 
This allows us immediately to write down the soft term renormalizations 

p'Xt = QYtA + YbAb + Y^^rriAs + 3a2mA2 + y^«i"^Ai, (92) 

16 7 
YtAt + &YbAb + YrAr + —aimA3+3a2mA2 + Yr^OiimAi, (93) 
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am 

(1) 
2 



I 1 



(1) 



^(1) 



-3^ 



5 

3 



Fi(m^ + ml + m^^ + A^) + Y,{m% + ml + m\^ + Al 



16 2 2 ^2 

-yagm^^ - ^a^m^^^ - — aim^^, 

16 



16 



m 



2 



o-t^ ^ 2 I 2 I 2 I /i2\ 16 2 4 

2y6(mQ + + m^^^ + A J - —a^mj,^ - — aim^^, 

2K(mi + m| + m^^ + A^) - yaim^^, 

3n(m^ + + m^^ + A^) + Yr{m\ + m| + m|^^ + A? 



which perfectly coincide with those of ref.(|T3||) 



(94) 

(95) 
(96) 
(97) 

(98) 
(99) 
(100) 
(101) 

(102) 
(103) 



6.2 Two and three- loop gaugino mass renormalization 

We calculate here the two and three loop gaugino mass renormalization out of a corre- 
sponding gauge /? functions. The RG /5 functions for the gauge couphngs in the MSSM 
are [l^ 



jk jf fg 



+ 



(104) 



where Yf means Yt, Y^ and Yj. and the coefficients 6j, Oj/-, bij^, aijf and aj/-„ are given 



in ref. 10 



For the gaugino masses we have 



^ hijkOijakimA, + + m^J - ^ aij/ajy/-(m^^ + m^^. - A/) 
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+ J2 (^ifgYfYgimA, - Aj - Ag 
fg 



+ 



(105) 



The same formulae can easily be obtained for the other soft terms. We do not write 
them down here due to the lack of space. Instead, the explicit formulae in the case 
when only top Yukawa coupling and are retained are presented in Appendix. 



Appendix. Three-loop renormalizations in the MSSM 

In this section, we present explicit formulae for rigid and soft term renormalizations 
in the MSSM in the three-loop approximation in the case when we retain only as and 
top Yukawa coupling Yf. 



The rigid renormalizations are [|10 



/?„3 = -3al + al{Ua3-4Yt) + al[^al-^a3Yt + 30Y,% (106) 

7, = (2Y,-^a,)-i8Y,' + ^al) + [i30 + 12Cs)Y,' (107) 

+ (f + 96C3)r^3 - (f + ^-fc3)YM + (^ + 320C3)ai], 

76 = --as--al + [--Ytal + i^ + 320C3)a% (108) 

7Q = (Yt - ^as) - i5Y,' + ^al) + [{15 + QC3)Y,' (109) 

+ (f + 48C3)y>3 - (y + ^^Cs)Y,al + + S20C,)a% 

jH, = (3Ft) - (9r,' - IGFtas) (110) 

+ [(57 + 18C3)F/ + (72 - 144C3)r/a3 " + ^(^C3)Ytal], 
16 . .__o 16 



r 16 16 

f3Y, = Ft|(6Fi-ya3)-(22r,2_i6F,a3 + _«2) +[(102 + 36C3)F,3 (111) 

+ ^F/a3 - (^ + 288C,)Y,al + (^ + 640C3)«^]} , 
/3^2 = /i^ {sFi - (9y/ - 16Fi«3) (112) 
+ [(57 + 18C3)y;' + (72 - 144C3)i;2a3 - (^ + 16C3)Ytal]^ , 

The corresponding soft term renormalizations read 

/9m3 = -3a3M3 + 28a^M3-4r,a3(M3- A) (113) 



104 



+ 347a^M3 - —aiYt{2M3 - A) + 30a3>^AM3 - 2A 



(3a, = (6FiA + ^^3^3) - [44r,% - 16yia3(A - M3) - ^«iM3] (114) 
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272 

+ [(306 + 108C3)y'M + ^r>3(2A-M3) 

- + 2^Kz)Ytal{A, - 2M3) - 3(^ + 640C3)«lM3], 

/3b = ZYtAt-[l^Y^At-l&Yta^{At-M^)\ (115) 
+ [(171 + 54C3)i;'A + (72 - 144C3)y/a3(2A - M3) 

- (l^ + 16C3ma^(A-2M3)], 

16 

13^2 = 2Yt{ml + ml + m]j^+Al)-—a^Ml-l&Y^\ml + ml + m]j^+2Al) 



16 

y 



a^Ml + 3(30 + 12C3)y;3(m2 + m| + m^^ + 3^2) 
16 



+ (y + 96C3)1^>3[(2A - M3)^ + 2{ml + m J + mf^J + M|] 

64 544 

- (y + — C3)i^t«^[(A - 2M,f + (m? + m| + m^J + 2M|] 
2720 

+ 12(— — + 320C3)Q;^M3^ (116) 



27 

16 , ,9 16 9 , ,9 
Pml = -—a^Ml-—alMl 



- -^yt<4Ut - 2M,f + {mi + m| + m^J + 2M|] 



3 

2720 

+ 12(^ + 320C3)aiM32, (117) 



16 

/3^| = rt(m2 + m| + m^^ + A2)-— a3M|-10F,2(m2 + m| + m|^ + 2A2) 
16 2,^2 



- y a^M3^ + 3(15 + 6C3)i;'(m,^ + m$ + m^, + 2,Ai) 
+ (y + 48C3)r/a3[(2A, - M^f + 2(m,2 + m| + m^J + M|] 

- (f + ^C3)l^t«^[(A - 2M3)^ + (m? + + m^J + 2M|] 



2720 

+ 12(yy + 320C3)«;M. (118) 

P^.^^ = Wt{ml + ml + m]j^ + Al)-lSY,\ml + ml + m]j^ + 2Al) 



+ 16ytQ;3[(A-M3)' + (m? + m| + m|J + M|] 

+ 3(57 +18C3)l^/(m,2 + mJ + m^^ + 3^2) 

+ (72 - 144C3)y/Q;3[(2^ - M^f + 2(m^ + m J + m^^J + M|] 

- (l^ + 16C3)rt«^p^_2M3)2 + (^2^^^ + ^2^J + 2M|], (119) 



where M3 = mA,- 
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Note added 

When this paper has already been finished we became aware of the paper where 
similar results were obtained. Our results coincide with their ones in many points. Ex- 
cept for the relation between the gauge coupling /5— function and the gaugino /3— function, 



the authors of Ref. []T6| have deduced the same relations between the renormalization 
group functions of the soft and rigid theories through the differential operators Di and 
D2 starting from the Yamada's rules 0. As for the Pa^PmA relation, their starting 



point was the Hisano-Shifman formula . Our approach is based on a consideration 
of the soft theory as a rigid one embedded into the external x— independent superfields, 
that are the charges and masses of the theory. The Yamada's rules together with the 
operator constructions Di and D2 are just the technical consequences of this approach. 
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